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ABSTRACT 

Powerful extended radio galaxies in the 3CR sample are observed out to 
redshifts of about 2. For redshifts greater than 0.3, the average lobe-lobe 
size of these sources decreases monotonically with redshift for all reasonable 
cosmological parameter choices. This suggests that the characteristic time for 
which an AGN produces highly collimated outflows that power radio emission is 
shorter for high-redshift sources than it is for low-redshift sources. The analysis 
presented here supports this conclusion. 

The relation between the active lifetime and the beam power of powerful 
extended radio galaxies is investigated here. It is found that the data are 
accurately described by a model in which the active lifetime of the source, t*, is 
written as a power-law in the energy extraction rate, Lj. The exponent of the 
power law is estimated to be (3 ~ 2.1 ± 0.6, where (3 is defined by t* cx LJ 13 ^ 3 . 
Note that the value of (3 for an Eddington limited system of zero is excluded 
by this analysis. The fact that (3 is constrained to lie within a certain range 
may be used to constrain models of large scale jet production and cosmological 
parameters. 

The comparison of the redshift evolution of characteristic source sizes with 
the average lobe-lobe size for powerful extended 3CR radio galaxies can be used 
to constrain cosmological parameters if three empirically estimated quantities 
can be accurately determined for a subset of the sources. As discussed here, one 
of these quantities, the lobe propagation velocity, is beset by potential biases 
that are not completely understood. The analysis presented here shows that 
these biases do not significantly affect the results on j3, but must be studied in 
more detail before cosmological parameters can be estimated precisely. Allowing 
for the potential biases mentioned above, best fits of the data yield a low of f2 D , 
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which is about la away from a flat, matter-dominated universe. 

Subject headings: cosmology:observations, galaxies:active theory — radio 
sources: galaxies 
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1. INTRODUCTION 

Active galactic nuclei (AGN) appear in a variety of forms which include powerful radio 
emitters produced by highly collimated outflows. The mechanism(s) that produces these 
outflows is(are) not known, though several models have been suggested (e.g. Blandford & 
Znajek f977; Blandford & Payne 1982; Wilson & Colbert 1995; Lightman & Eardley 1974; 
Rees et al. 1982; Narayan & Yi 1994, 1995). 

Powerful extended radio galaxies can be used to probe AGN central engines, as 
described here and by Daly (1994, 1995). Radio sources are referred to as "powerful" 
if their radio powers at 178 MHz satisfy Pn$ > 3h~ 2 x 10 26 W Hz -1 sr _1 , assuming a 
deceleration parameter q Q = and parameterizing Hubble's constant in the usual way: 
H = 100 h km s" 1 kpc" 1 . These sources are all on the FRII (edge-brightened) side of the 
FRI-FRII break defined by Fanaroff & Riley (1974). They are "extended" because their 
core-lobe separations range from about 25 h~ x kpc to 240 Yr x kpc (e.g. the radio galaxies 
listed by Wellman, Daly, & Wan 1997a and 1997b; hereafter WDW97a and WDW97b), 
which indicate that the outflows are interacting with intergalactic or intracluster gases. 
Only galaxies, and no quasars, are considered here to minimize projection effects, since it is 
widely believed that powerful extended radio galaxies have lobes and bridges that lie closer 
to the plane of the sky than those of quasars (see §3.1 below, and Wan & Daly 1997a). 

One key to the properties of collimated outflows is the relation between the rate at 
which energy is channeled into the outflow, known as the beam power or luminosity in 
directed kinetic energy, and other parameters such as the total time the outflow exists, 
and the energy available to power the outflow. In an Eddington limited system, the rate 
of energy extraction, Le, is proportional to the mass of the central compact object, M. 
Since the energy available to power the source is equal to the mass of the central compact 
object times the emission/accretion efficiency, E ~ i]M, the total time the outflow exists is 
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ts — E/L E oc rjM/M oc T]. Thus, the total time the outflow exists for an Eddington limited 
system depends only on rj, and not explicitly on L E or E. 

The relation between the beam power and the total time an AGN produces collimated 
outflows is empirically investigated here for the case of powerful extended radio galaxies. 
The data can be understood if the total time the AGN produces collimated outflows 
depends on the beam power, and thus, the outflows are not Eddington limited systems. 
The empirical relations obtained here can be used to constrain models of energy extraction 
from AGN, and may provide insight on the conditions at the very core of AGN. 

The method and model described here can be used to constrain cosmological parameters 
if the redshift behavior of empirically determined quantities can be estimated to relatively 
high accuracy (Daly 1994). The constraints on cosmological parameters that can be placed 
with the existing data are described in §||. Modest constraints (i.e. about 2 a) on a flat, 
matter-dominated universe may be placed at the present time. More precise constaints 
can be placed when the lobe propagation velocity is determined more accurately, or can 
be observed directly. As described in §4.3| , uncertainties that affect the lobe propagation 
velocity do not affect, within the errors, the estimate of the model parameter /3, which 
relates the beam power and the total time an AGN produces collimated outflows. 

A basic model that describes the observed characteristics of powerful extended radio 
sources rather well is presented in §2. The radio source samples used in the present study, 
including their limitations, are described in §3. The application of the data to the model, 
and constraints on the relation between the beam power and the time for which the outflow 
occurs are detailed in §4. The use of current data to estimate and constrain cosmological 
parameters is discussed in §5. A general discussion follows in §6. 
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2. A MODEL OF THE LOBE SEPARATION IN POWERFUL EXTENDED 

RADIO SOURCES 

2.1. Foundations of the Model 

The foundation of a model which describes the lobe separation of powerful extended 
radio sources is discussed in detail by Daly (1994, 1995). Powerful extended radio sources 
have regular, straight bridges and lobes that are thought to propagate supersonically 
outward into the ambient gas surrounding them ([Leahy fc Williams 1984J , hereafter LW84; 
Alexander fc Leahy 19871 , hereafter AL87). The lobe propagation velocities are computed on 



the basis of spectral aging along the bridges (e.g. Myers & Spangler 1985; Leahy, Muxlow, 
& Stephens 1989, hereafter LMS89; WDW97a,b). 

It should be noted that results from standard spectral aging have been criticized 
recently (Katz-Stone, Rudnick, & Anderson 1993; Rudnick, Katz-Stone, & Anderson 1994; 
Rudnick & Katz-Stone 1996; Eilek & Arendt 1996; Eilek 1996). These authors have argued 
that syncrhotron aging may not be the correct or only explanation for the spectral curvature 
observed in the bridges of extended radio sources. Until more conclusive data are obtained 
to settle this controversy, we adopt the practice of using lobe velocities from spectral 
aging analysis since WDW97b find that these velocities and the inferred surrounding gas 
temperatures agree well with results from lobe assymetry analysis and X-ray observations. 

WDW97a find the data agree well with a model in which the bridges of these sources 
expand laterally and the plasma inside undergoes adiabatic expansion. Their analysis 
suggests that the backflow velocity of the plasma within the bridge is small compared to 
the lobe advance velocity and can be neglected. It should be noted that backflow can 
be negligible along quiescent bridges, but significant near the hotspot, in agreement with 
numerical simulations (e.g. Norman 1996). 
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The ambient gas density around the lobes of powerful extended radio sources can be 
estimated on the basis of ram pressure confinement since the lobes propagate supersonically 
(see for example Daly 1994, 1995; WDW97b). It can be shown that the number density of 
the ambient gas, n a , is given by 

n a oc f(b)3g* , (1) 

where vl is the lobe propagation velocity, B m i n is the minimum energy magnetic field, h is 
the offset of the lobe magnetic field, Bl, from minimum energy magnetic field (Bl = bB m i n ), 
and f(b) = (4/3)6~ 3 / 2 + b 2 . WDW97a compute values of b for Cygnus A by comparing 
their estimates of pressure in the lobes, pressure in the bridge, and mach number of lobe 
advance to estimates based on X-ray observations of gas density, thermal pressure and 
temperature of the surrounding ambient gas. In particular, each of these three independent 
estimates gives b ~ 0.25 for Cygnus A (see §3.2. lp , in agreement with the result of Carilli et 
al. (1991). WDW97b show that the source to source dispersion in b must be rather small, 
less than about 15%. 

Daly (1990) derives an expression for the beam power or the luminosity in directed 
kinetic energy, Lj, by equating the work done by the lobe while it propagates into the 
ambient gas to the energy supplied by the collimated outflow: 

Lj oc n a a\vl oc (v L /kf, where k = (n a a 2 L Y 1,i (2) 

and cll is the lobe width. It should be noted that the typical Lj for radio galaxies in the 
sample examined below (§ |3.2| ) is about 10 45 ergs s _1 assuming b = 0.25, and it appears 
that Lj is time independent for a given source (Wan, Daly, & Wellman 1996; Wan & Daly 
1997c). For a roughly constant Lj, the energy supplied by the central engine to power 
the collimated outflows is Ei ~ Lj t* where i* is the lifetime of the collimated outflows. 
For t* between 10 7 to 10 8 , which agrees with spectral aging (e.g. AL87; LMS89; Liu, 
Pooley, and Riley 1992, hereafter LPR92), the energy supplied by the central engine is 
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Ei ~ 10 60 ergs ~ lO 6 M c 2 . For emission/accretion efficiencies of 0.01 to 0.1, the mass of 
the central compact object would be 1O 7 M to 1O 8 M which is similar to those discussed 
in models of jet formation in powerful FRIIs (e.g. Wilson & Colbert 1995). 



2.2. Relation Between Active Lifetime and Beam Power 

The data clearly show that the average lobe-lobe size of powerful extended 3CR 
(Bennett 1962) radio galaxies decreases with redshift for z > 0.3, while the lobe propagation 
velocity computed on the basis from spectral aging tends to increase with redshift and radio 
power (see Figure [l| in this paper and Figure 19 in WDW97a). Also, it appears that there is 
no relation between radio power and lobe-lobe size at fixed redshift (e.g. Lacy et al. 1993; 
Nesser et al. 1995; Wan & Daly 1997b). Since the lobe-lobe size is proportional to the lobe 
propagation velocity times the lifetime, it can be inferred that high redshift 3CR galaxies 
have shorter lifetimes than low redshift sources. The higher redshift sources also have larger 
beam powers than the lower redshift sources since radio power and beam power are roughly 
proportional (Wan, Daly, & Wellman 1996; Wan & Daly 1997b). A shorter lifetime for a 
more powerful source is contrary to what is expected in an Eddington limited outflow, since 
the Eddington lifetime depends only upon certain efficiency factors and is independent of 
the beam power and the total energy of the central compact object. Though it should be 
noted that such an outflow could still be Eddington limited if the efficiency factors vary 
with the beam power. 

Let the total time that an AGN produces two highly collimated oppositely directed 
outflows be t*, and let the rate of growth of the radio bridges, vl, be roughly constant over 
the source lifetime, which is supported by spectral aging analyses (AL87, LPR92, Daly 
1994, 1995). The average or characteristic lobe-lobe size a source would have if it could be 
observed over its entire lifetime is ~ vl t*. Following Daly (1994), £* is related to the 
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energy extraction rate or beam power, Lj, by a power law: 

UocLj-W . (3) 

This implies that the characteristic source size is 

D, oc k Lf~ P)l * oc v\~ p , (4) 

where k = (nad^)" 1 / 3 (see eq. Clearly, for (5 — 1, the average source size depends 
only upon n a and ol, and is independent of Lj, whereas for /3 > 1, the average source size 
decreases as the beam power and radio power increase. 



The implications of eq. (|3]) are far-reaching. It implies that the total energy available 

j t* oc Lj 



to power the outflow is fixed at some initial value, since Ei ~ Lj t* oc L^ /3 ^ 3 ' ) , assuming 



that Lj is roughly constant over the source lifetime which is supported by the current data 
(see §2.1). Note that for (3 ~ 3, the energy extraction rate Lj is independent of the energy 
Ei available to power the outflow, while for (3 ~ 0, the beam power and AGN lifetime are 
independent and Lj oc Ei, analogous to an Eddington luminosity. 

In practice, D* is estimated using the expression introduced by Daly (1994), originally 
referred to as Z*: 

/ 1 \ 2/3/3 

D* oc ( -J— ) . (5) 



The normalization is chosen so that at 2 ~ matches the observed average lobe-lobe 



size for sources at this redshift that are subject to the power cut discussed in § [3.1| . The 
magnetic field strength, Bl, in a synchrotron-emitting radio source and the lobe propagation 
velocity based on spectral aging (Myers & Spangler 1985), vl, can be estimated using the 
field strength that minimizes the total energy in relativistic electrons and magnetic field 
(e.g. Miley 1980). The estimates of Bl and vl can also be computed using magnetic 
fields systematically offset from the minimum energy field, which is discussed in detail by 



WDW97a and WDW97b (see pXID 
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Assuming the energy density of the field in the radio bridge is large compared to 
that in the microwave background radiation (i.e. synchrotron losses dominate over inverse 
Compton losses) and a spectral index of a = 1 (S u oc z/~ a ), it can be shown that 

D, oc {a rf A ' 7 ' 2 ^ (1 + z )^~^) , (6) 

where z is the redshift of the radio source, and (a a r) is the coordinate distance to the 
source. The value of (3 indicated by the data for all parameter choices is (3 ~ 2 (see j j4.3|) , 
and for this value, D* oc (a r) -0 - 8 . The ratio < D > /£?*, where < D > is the average 
lobe-lobe size, has a rather strong dependence on cosmological parameters via (a Q r): 
< D > jD m oc (a D r) L8 (1 + z)' 1 , for (3 ~ 2. The results presented in subsequent sections 
include the effects of inverse Compton scattering with microwave background photons and 
the observed spectral indices, but the assumptions made above show the dependence of 
on cosmology in one approximation. 

The main hypothesis of the model is that single epoch radio data can be used to 
estimate the average lobe-lobe size the source would have if it could be observed over its 
entire lifetime, called D*. Since very powerful FRII radio galaxies form a very homogeneous 
population, the average size of a given source should be close to the average size of the 



population at the same redshift, as discussed in detail in § [4.2| . Thus, one test of the model 
is whether tracks < D > independent of redshift for a large sample of sources selected 
using similar criteria; this is shown to be the case in §4.3| . If the sources are sampled 
randomly during their lifetimes, the distribution of t/t* should be redshift independent, 
where t is the age of the source when it is observed and t* is the lifetime of collimated 



outflows; this is shown to be the case in §4.2|. Also, if the model correctly describes the 



data, then the ratio < D > / D* should exhibit much less scatter than the ratio D/D*; this 
is shown to be the case in 
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3. RADIO GALAXIES EXAMINED 
3.1. The Sample Used to Estimate < D > 



In order to able to use the model described in §2^, the sources examined must have 
lobes that propagate supersonically, and have no significant backflows or bridge distortions. 
Thus, only FRII radio galaxies with Pn$(q — 0) > 3h~ 2 x 10 26 W Hz -1 sr -1 are examined 



here, (see §1 & |2.1|) . Radio-loud quasars are not included in this work since the lobe-lobe 



sizes of radio-loud quasars evolve differently with redshift than radio galaxies (e.g. |Singal 



1988|) , their bridges are more distorted than those of radio galaxies (LMS89), and they may 



suffer from more serious projection effects. 

The 3CR radio galaxies with galactic latitude greater than 10° are completely identified 
and all have spectroscopic redshifts. On the basis of published 178 MHz fluxes and spectral 
indices compiled by Spinrad et al. (1985), 81 radio galaxies from the list compiled by 
McCarthy, van Breugel, & Kapahi (1991) have Pi 78 (g = 0) > 3/i" 2 x 10 26 W Hz^sr" 1 . 
These sources are used here to estimate < D >. Cygnus A (3C405), which is not listed by 
McCarthy et. al. (1991) due to its low galactic latitude (6°), is added to this sample since 
it is the closest high power FRII radio galaxy (z = 0.056), is well studied (e.g. Carilli et al. 
1991), and is included in the sample described below (§ p.2|) . 



Table [| gives the mean and median lobe-lobe sizes for the high-power FRII 3CR 
galaxies assuming a flat, matter-dominated universe with (Q D = 1.0, Qa — 0), an open, 
curvature-dominated universe with (Q = 0.1, fi a = 0), and a spatially flat universe with a 
nonzero cosmological constant (fl = 0.1, Oa = 0.9); here and throughout Q Q refers to the 
ratio of the mass density to the critical value at the current epoch, refers to the ratio 
of the energy density of a cosmological constant to the critical density, = 1 — fi G — 
is the normalized curvature. A value of h — 1 is adopted when a value must be chosen, 
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though h scales out of many quantities, and the dependence on h is typically weak when it 
does enter. 

Figure |l| shows the mean lobe- lobe sizes from Table 1. Excluding the lowest redshift 
bin, an obvious decrease in the physical size with increasing redshift can be seen for the 
cosmological parameters considered. Note that the lowest redshift, which contains Cygnus 
A, is not an obvious outlier in the comparison of and < D > as discussed in §4.3 below 
and by Guerra & Daly (1996). 



3.2. The Sample Used To Estimate 

An estimate of requires estimates of the lobe radius, the lobe magnetic field 
strength, and the lobe propagation velocity of the radio source (see eq. |5]). Fourteen 3CR 
radio galaxies have enough radio bridge data published to estimate these parameters, 



and satisfy the power cut described in § [2.1| . Wellman and Daly (1996), and WDW97a,b 
reanalyzed the radio maps from LMS89 and LPR92 for their studies of the bridge structure 
and gaseous environments of powerful extended radio sources, and computed vl, Bl, and cll 
for both data sets in a similar manner. Table 2 lists the 14 radio galaxies from WDW97a,b 
all of which are used here. Note that since they are all high-power 3CR radio galaxies, they 
are all included in the larger, full high-power 3CR sample of radio galaxies, described in 



§ |3TT| . One bridge from 3C239 was excluded by WDW97a,b because its morphology suggests 
sideflow (see WDW97a,b for details). WDW97b compute values for the lobe half-width, 
a,L, measured 10 h^ 1 kpc from the hotspot toward the host galaxy, the magnetic fields 
measured 10 h~ x kpc and 25 h~ x kpc from the hotspot toward the host galaxy, and the lobe 
propagation velocity, vl, which is affected by offsets in the magnetic field from minimum 
energy conditions and corrections for redshift evolution of spectral index as discussed below. 
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3.2.1. Uncertainties that Affect Estimates of 

The magnetic field in the radio lobe, Bl, enters into the estimate of the lobe pressure, 
Pl, and the magnetic field strength of the radio bridge, E>b, and Bl enter into the estimate 
of the lobe propagation velocity. Both values affect estimate of the ambient gas density, 
n a (WDW97b). Several observations indicate that the magnetic field strength in powerful 
extended radio sources is less than the minimum energy magnetic field. Expressing the 
true magnetic field as B = bB min , where B min is the minimum energy magnetic field, 
Carilli et al. (1991) find that b = 0.3 in order for ram pressure confinement of the lobes 
of Cygnus A to be consistent with X-ray measurements of the ambient gas density. Perley 
Sz Taylor (1991) find a similar value for 3C295 based on ram pressure confinement of the 
lobes, while Feigelson et al. (1995) and Kaneda et al. (1995) find similar values for Fornax 
A by comparing the radio emission with X-ray produced by inverse Compton scattering 
of microwave background photons with the relativistic electrons that produce the radio 
emission. WDW97a use 3 independent and complementary methods to estimate b in 
Cygnus A, and their results agree with the Carilli et al. (1991) result. WDW97b also show 
that the source to source dispersion in b must be small, less than about 15%. However, 
both b = 0.25 and b = 1.0 cases are examined below for completeness (§[4.3|, Tables 3-5). 

Another uncertainty that affects is whether the evolution of the radio spectral 
index with redshift introduces a systematic error on the lobe propagation velocity. It has 
been noted for some time that the radio spectral index of the 3CR sample and of other 
samples increases systematically with redshift; this could be due to spectral curvature or 
due to other causes (e.g. Rottgering et al. 1994). The radio spectral index is an important 
ingredient in estimating the lobe propagation velocity via the effects of spectral aging 
of relativistic electrons. Though it should be noted that for (3 ~ 2, which is indicated 
by the data independent of systematic effects on the lobe propagation velocity (jfPf), 
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D* oc (ol-Bl) -4 / 3 v l , which is a fairly weak dependence on vl- It is not clear whether 
the data should be corrected for the systematic increase of the radio spectral index with 
redshift. This correction, referred to as the a-z correction, does not change the low-redshift 
velocities, but decreases the high- redshift velocities by, at most, a factor of two (see 
WDW97a, Figure 19). Results obtained both with and without the a-z correction are 
examined below (jp~3|, Tables 3-5). 

4. EXAMINING THE MODEL 

4.1. Computed -D*'s 

Each bridge has a characteristic core-lobe length, r* ~ D*/2, where is obtained 
using eq.(|5|) and inputing the values of vl, o>l, and Bl from each bridge. Note that one 
can estimate using the average of the input parameters in eq. (|5]) over both bridges, 
or by adding r* for both bridges. Both methods would be equivalent if the sources were 
completely symmetric, but the latter method accounts for any slight asymmetry in these 
sources and is used here. As an example for one set of parameter choices, Table |2| lists r* 
computed using the results for the radio galaxies from WDW97b, and assuming b = 0.25 
without the a-z correction, (Q = 0.1, ^a = 0), and (3 = 2.0. This example is chosen because 
(Q = 0.1, = 0) is the moderate choice of the three example cosmologies introduced in 
Table 1, and (3 = 2.0 is the value of (3 indicated by the data for all parameter choices ( §4.3|) . 
Estimates of throughout this paper are computed by adding both r* values, except for 
3C239 which has r* for only one bridge and has -D* set equal to twice the available r*. The 
normalization of is chosen so that the ratio of the lobe-lobe D* for Cygnus A (3C405) 
to < D > for the lowest redshift bin is unity. 

Figure ^|a shows as a function of (1 + z) assuming (Q Q = 0.1, ^a = 0), b = 0.25 
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without the a- z correction, and (3 = 0.0 which is indicative of an Eddington limited system; 
Figure 0b shows the same except it is assumed that (3 = 2.0 which is the value indicated 



by the data for all parameter choices in §4.3. In these figures and all other figures in this 



work, diamond symbols represent radio galaxies from the LMS89 sample, and star symbols 
represent radio galaxies from the LPR92 sample. These figures show that estimates 
for LMS89 and LPR92 sources agree well with each other for the range of (3 considered, 
although the 2 samples are more similar for (3 = 2 than for (3 = 0. Note that the LPR92 
sources have smaller angular extent than the LMS89 sources (Table 2), and have smaller 
lobe-lobe sizes which implies they are younger (see § [4.2| below). The fact that the smaller 
sources from LPR92 and the larger sources from LMS89 give consistent at similar 



redshifts, indicates that the model discussed in §2.2 is yielding a characteristic length that 



is independent of when the source is observed during its lifetime. 

Figure 0a shows that increases as redshift increases for (3 = 0.0, which differs from 
the redshift dependence of the average lobe- lobe sizes (see Figure [l]). Figure ||b shows 
that, ignoring the lowest redshift bin (Cygnus A), D* decreases with increasing redshift for 
(3 = 2.0. The general trend is that increasing (3 makes the slope of vs. (1 + z) decrease 
(i.e. become more negative) excluding the lowest redshift bin. Thus, for (3 = 2.0, the 
redshift behavior of tracks that of < D >. 

For (3 = 0.0, -D* for 3C427.1 (z = 0.572) is similar to the for Cygnus A (z = 0.056), 
while for (3 = 2.0, .D* for 3C427.1 is similar to those for the rest of the 0.3 < z < 0.6 sources. 



Excluding this source slightly decreases, by about 0.3, estimates of (3 (see § [4.3| ). Cygnus 
A (z = 0.056) has low values for when compared with most sources in the redshift 
interval 0.3 < z < 0.9. This is intriguing since < D > has a relatively low value in the 
lowest redshift bin and tracks the redshift behavior of when Cygnus A is included (see 
Figure [l]). Excluding Cygnus A does not affect the estimates of (3, but could slightly affect 
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estimates of Q a using this data set (see Guerra & Daly 1996). 

4.2. Comparing D to D* 

The observed lobe-lobe size of a powerful extended radio source is D ~ 2vr,t where t 
is the age of the source when observed, and vl is the lobe propagation velocity which is 
assumed to be roughly constant over the lifetime of a source. An alternate way of expressing 
the size for a given source is 

D = 2Q.)D.. (7) 

In a given redshift range (z — 5z to z + 8z) of a sample of radio sources, the distribution of 
.D* values should be peaked around some central value D* z in a given redshift range (z — 5z 
to z + 5z). Taking the values of in this redshift range to be equal to the constant D* z , 
the mean lobe-lobe size is 

( D )= 2 ((£) D -H(^) c ~ (8) 

If we assume for a given sample that the distribution of t/t*, the fraction of a source's 
lifetime at which it is observed, is a constant over £/i*, then < D >~ D^. More generally, 
if the t/t* distribution of a sample does not depend on redshift, then < D >oc D* z , and the 
ratio < D >/D* z should be constant and independent of redshift. It is this premise that 
allows an estimate of (3 and cosmological parameters. 

Figure |3] shows the ratio of the physical size D to the characteristic size for the 
sample of radio galaxies examined here, assuming b = 0.25 and not including the a-z 
correction. In Figure [| (Q a = 0.1, = 0.0) is assumed since it is the moderate choice of the 
three example cosmologies in Table 1, and (3 = 2.0 is assumed since it is the value indicated 
by the data for all other parameter choices (see |Q| ). Since eq. (0) gives t/t* oc D/D* for 
a given source with constant vl, the sample shows no significant t/t* redshift evolution. 
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Fitting D/D* = a(l + z) n we obtain a = 1.15 ± 0.07 and n = —0.04 ± 0.09; similar results 
are found for b = 1.0 and including the a-z correction. Thus, the data are consistent with 
the sources being observed at random times during their lifetimes. The reduced y 2 of the 
fit is 41 which is due to the large scatter of D/D*. It should be noted that the scatter in 
D/D* is much greater than that of < D > / D* (see Figure |8|b), which indicates that though 
the sources are sampled over a broad range of fractional ages, sources at a given redshift 
yield similar D*. 



4.3. Constraints Obtained by Comparing < D > to at a Given Redshift 

A quantitative constraint on the parameter j3 may be obtained by fitting the ratio 
< D > /D* to a constant, independent of redshift, and finding the value of f3 that minimizes 
the reduced x 2 ■ Table || lists the values of (3 that minimize the reduced x 2 ■, Pmin, for 
(Q = 1.0, Q\ = 0), (Q = 0.1, = 0), and (Q = 0.1, = 0.9). Four cases are examined 
in Table (A) b = 0.25 without the a-z correction, (B) b = 0.25 with the a-z correction, 
(C) b = 1.0 without the a-z correction, and (D) b = 1.0 with the a-z correction. Error 
estimates for the fitted parameters here and thoughout are the 68% confidence intervals 
that correspond to the regions of the parameter space within which the \ 2 increases by no 
more than 1.0 from the minimum value (e.g. Press, Teukolsky, Vettering, & Flannery 1992; 
§15.6). Although the errors that enter into computing the x 2 ma y n °t be Gaussian, we shall 
assume Gaussian errors as a rough approximation in order to use the x 2 as an estimate of 
errors for fitted parameters. The /3 m i n 's from Table ||] are all consistent within 2a for all the 
cases (A-D) and the three cosmologies examined. A range of about 1.5 to 2.6 is found for 

min- 

Figures £|a &j|b show contour plots of the reduced x 2 as a function of (3 and Q a , for 
case A described above, assuming ^a = and Qk = respectively. It is clear from these 
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figures that the constraints on (3 are not significantly affected by the assumed cosmology, 
and a consistent range of f3 min from 1.5 to 2.75 emerges. Most significantly, (3 = gives a 
reduced \ 2 just slightly greater than 6, for all cosmologies shown. Thus, taken at face value, 
the data are not consistent with Eddington limited outflows in powerful extended radio 
sources. Contour plots for the cases B, C, and D described above are also shown in Figures 
[5]a-b, |5]a-b, and [j]a-b respectively. The constraints on (3 are not significantly different for the 
four cases (A-D) examined, and all cases give a high reduced x 2 , greater than 3, for (3 = 0. 
Note that these figures only depict the range of Q a from to 1, but fits have been extended 
outside this range, including Q Q < which gives divergent values of the coordinate distance, 
(a Q r), for k = (see Table (|). 

Table ^ shows the best fit values of Q D and (3 for a constant < D >/D* independent 
of redshift. The fits for Q vary by about 0.2 to 0.4 depending on the case (A-D) described 
above, which follows from Figures A value of Q Q less than unity is indicated by these 
fits for all parameter choices, and Q = 1 is inconsistent at about the 2a level. These 
constraints on Q Q are as strong as those from studies of large scale velocity fields which 
indicate that fl a > 0.3 at about the 2Aa level (e.g. Dekel et. al 1993, Dekel & Rees 
1994, Hudson et. al 1995). The use of radio powerful extended radio galaxies to constrain 
cosmological parameters is discussed further in §|5| below. 

Table |^ shows the results of fitting < D to a power law in (1 + z), 

<D>/D*(x(l + z) p . A value of p = is expected for the correct parameter choices (see 
§ |4.2j ). Best fit values for (3 and p are shown for the four cases described above (A-D) and 
three choices of cosmological parameters. It is interesting to note that the fits for (3 converge 
to similar values for different cosmologies when a power law fit is allowed. The fits where 
(Q = 1.0, Q\ = 0) are the only ones that give a significantly non-zero value of p, in which 
case p is about 2a less than zero. Figures |||a, |||b, and |c show the ratio < D >/D* and the 
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best fit to a constant (the weighted mean in log-space), assuming b = 0.25 without the a-z 
correction and (3 = 2.0, for (0 = 1.0, ft A = 0), (Q Q = 0.1, fi A = 0), and (fi = 0.1, fi A = 0.9) 
respectively. The slopes found in the power law fits for (Q a = 1.0, Q\ = 0) are apparent 
in Figure [8]a, and can be explained by the incorrect choice of cosmological parameters. 
This implies that if the universe is flat and matter-dominated, then there is some redshift 
evolution that has not been accounted for. 

The possibility that certain sources in the sample may not be representative of 
powerful extended radio sources at their respective redshifts was considered. For example, 
the only low redshift source for which we have an estimate of is Cygnus A (3C405, 
z = 0.056); perhaps Cygnus A is an unusual source. A source that is similar to Cygnus A in 
many respects, including size, lobe propagation velocity, surrounding ambient gas density, 
and ambient gas temperature, is 3C427.1 (WDW97a,b). Also, two sources, 3C68.2 and 
3C239, do not seem to follow the predictions of adiabatic bridge expansion, according to 
WDW97a. All of the analyses described above have been repeated while excluding different 
combinations of these sources (3C405, 3C427.1, 3C68.2, and 3C239), and the differences 
are negligible except for a slight reduction of (3 min by about 0.3 when 3C427.1 is excluded 
holding all other parameters fixed. The reduction of f3 m i n when 3C427.1 is excluded occurs 
because lower (3 tends to make 3C427.1 an outlier (see Figures |2|a,b). 

The data give an estimate of (3 that is insensitive to virtually all other parameter 
choices and data cuts. A reasonable estimate of (3 is (3 = 2.1 ± 0.6 based on Figures |]-[7| and 
Tables 3-5. It is safe to state that (3 must be between 1 and 3, and that these outflows are 
not Eddington limited, (3 ^ 0. The main conclusion that can be drawn about cosmological 
parameters is that a flat, matter-dominated universe (0 G = 1.0, fi A = 0) is about 2a away 
from the best fit values (see §|5| below). 
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5. CAN POWERFUL EXTENDED RADIO SOURCES BE USED FOR 

COSMOLOGY? 



In order to be able to use the model outlined in § [2.2| to estimate cosmological 
parameters, accurate estimates of or,, Bl, and Vl are needed (see eq. Though two 
inputs to Vi are not completely understood at this time, it is still interesting to apply the 
method to this data set and consider all possible effects on vl to see the implications for 
cosmological parameters. 

One factor that is not completely understood is the offset of magnetic field strength 
from the minimum energy value. This enters as a scaling factor of Bl and only affects 
the redshift behavior of D* through the redshift behavior of vl- It turns out that the 
dependence of D* on vl is rather weak for (3 ~ 2 (see eq. RD: oc (a^i^) -4 / 3 t>]/ 3 , but the 
effect of the offset from minimum energy conditions on vl is considered nonetheless. As 
the magnetic field strength decreases from the minimum energy value, the role of inverse 
Compton cooling begins to become important relative to the role of synchrotron cooling 
in these sources. The affects of both types of cooling on the radio spectral index is used 
to estimate spectral ages, and thus vl- As discussed by WDW97b, for b = 0.25, the 2 
cooling mechanisms are of comparable importance for many sources at high redshift. This 
offset affects both the redshift behavior of the vl (shown in Figure 4 of WDW97b), and the 
dependence of on (a D r), the coordinate distance to the source. 

The second effect that may change estimates of the lobe propagation velocity is the 
increase of the radio spectral index with redshift. It is not clear whether the data should 
be corrected to account for the observed systematic increase of the radio spectral index 
with redshift (discussed in detail by WDW97b). For 3CR radio galaxies, vl increases with 
redshift, as do the radio power and the radio spectral index. This increase is well known 
and has been noted by many authors. If, for example, the initial radio spectrum is not a 
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power law but has some curvature (see e.g. Rottgering et al., 1994), the observed spectral 
index should be corrected for the observed redshift evolution as was done for the a corrected 
data (WDW97b). On the other hand, if the correlation is not due to, or only partly due to, 
spectral curvature than such a correction is not appropriate. 

To account for these two uncertainties, both b = 0.25 and 6=1 are considered, and 
both ol-z corrected and a-z uncorrected velocities are considered (see § |4.3|) . The results are 
not strongly dependent on these changes (see Table £|), as expected since -D* oc v^ 3 . The 
effect of using a value of b = 1.0 instead of b = 0.25 is to increase the best fit Q Q by about 
0.2 and the effect of including the corrections to spectral index is to lower the best fit Q Q 
by about 0.2 (see Table A precise fit for Q Q cannot be made until the parameters that 
are varied between the four cases (A-D) described in §4.3 are known. However, the current 
range of Q Q allowed can be discussed. 

An ideal solution to the uncertainties described above would be to identify an 
independent estimate of the lobe propagation velocity for all or some of these sources. This 
could be compared with the lobe propagation velocity estimated using the spectral aging 
model used here. The velocity can be independently estimated if either the ambient gas 
temperature or the ambient gas density can be estimated. WDW97a have shown that 
the geometrical shape of the radio bridge may be used to estimate the Mach number of 
the lobe, so the combination of the Mach number and the ambient gas temperature can 
be used to estimate the lobe propagation velocity. Perhaps observations by AXAF will 
be used to measure the ambient gas temperature since its spatial resolution and spectral 
coverage could allow it to separate X-rays produced by the AGN from those produced 
by thermal bremsstrahlung from the hot ambient gas. Daly (1994, 1995) and WDW97b 
have shown that the lobe propagation velocity and lobe pressure may be used to estimate 
the ambient gas density; if the ambient gas density can be estimated independently, the 
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lobe pressure and ambient gas density may be combined to solve for the lobe propagation 
velocity. Studies of the gaseous environments of these sources using the Sunyaev-Zeldovich 
effect or high resolution AXAF measurements may lead to independent estimates of the 
ambient gas density, which may be combined with the lobe pressure to solve for the lobe 
propagation velocity. 



6. DISCUSSION 

A relation between the active lifetime during which an AGN produces a highly 
collimated outflows, t*, and the beam power Lj, is introduced which can reconcile the 
observed evolution of radio source size with redshift, and the independence of radio source 
size and radio power at a given redshift for powerful extended radio sources (see §2.2). 



What emerges from this model is a characteristic length scale, D*, that can be thought of as 
the lobe-lobe size averaged over a source's lifetime, which can be predicted from measurable 
quantities such as the lobe width, the lobe propagation velocity, and the lobe magnetic field 
strength. 

The model parameter /3, which relates t* and Lj, is estimated by comparing the redshift 
evolution of to that of the average lobe-lobe sizes for powerful extended radio galaxies. 
For all the cases and cosmologies considered, a consistent range of (3 emerges (§ [4.3| ). The 
best fit values allow (3 to range from about 1.5 to 2.75, so a good rough estimate of j3 is 
2.1 ± 0.6. This is consistent with the range of (3 = 1.5 ± 0.5 obtained by Daly (1994, 1995). 

These data are not consistent with (3 = 0, which is expected for an Eddington limited 
system (t* independent of Lj as discussed in §|2.2|). In Figure 0a, where (3 = 0.0 is assumed, 
the estimates are clearly increasing with redshift, which is opposite the evolution of 



< D > for the high-power 3CR sample described in §3.1. The data are best fit with 
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(3 > 1.0 which yield reduced x 2 ' s l ess than 1.5, whereas fits assuming (5 = yield reduced 
% 2 's greater than 3 (see Figures §-0). The values of /? obtained here imply that .D* should 
decrease for increasing Lj or vl (see § |2.2|) , which is consistent with the redshift evolution 
of the average sizes of high-power 3CR sources (Figure |l|), and the redshift evolution of Lj 
and vl (discussed by Wan & Daly 1997b). 

A flat, matter-dominated universe (Q Q = 1.0, Qa = 0) is allowed at about the 2a 
level for the full range of parameter allowed. It should be noted that the significance 
of the constraints on Q obtained here are similar to those obtained from the study of 
the large scale velocity fields (e.g. Dekel et. al 1993, Dekel & Rees 1994, Hudson et. 
al 1995). Stronger cosmological constraints can be placed using this method when the 
lobe propagation velocity can be more accurately determined, as discussed in §|5|, or can 
be measured directly by comparing X-ray temperature measurements with the shape 
of the radio bridge, as described by WDW97a, or by comparing the X-ray density or 
Sunyaev-Zeldovich pressure with the pressure of the radio lobe, as described by WDW97b. 
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Table 1. Lobe-Lobe Sizes for P 178 > 3h~ 2 x 10 26 W Hz" 1 sr" 1 FRII 3CR Galaxies. 









n Q = i.o, ^ A = o.o 


n Q = 0.1, ^ A = o.o 


n = 0.1, ^ A = 0.9 






no. 


< D > 




< D > 




< D > 




bin 


z range 


sources 


(kpc/h) 


(kpc/fc) 


(kpc/fc) 


(kpc/h) 


(kpc/fc) 


(kpc/h) 


1 


0.0-0.3 


3 


66±14 


52±14 


68±14 


55±14 


72±13 


60±13 


2 


0.3-0.6 


15 


177±42 


160±43 


197±47 


175±48 


227±54 


200±55 


3 


0.6-0.9 


27 


126±18 


127±18 


149±21 


147±21 


179±25 


176±26 


4 


0.9-1.2 


19 


91±21 


51±16 


114±27 


63±21 


141±33 


79±26 


5 


1.2-1.6 


10 


83±30 


39±18 


110±40 


51±25 


137±50 


63±31 


6 


1.6-2.0 


8 


50±17 


37±18 


71±24 


53±25 


88±30 


66±31 
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Table 2. Radio Galaxies with Presented Here. 









Log[Pi 78 (g = 0)J 


n 
V 


Map 








c 

source 


z 


bin 


/ 7 — 2 TT — 1 — 1\ 

[n W Hz sr J 


( \ 

(arcsecj 


t ) r 

Kel. 


(kpc/ft) 


(kpc/ft) 


(kpc/ft) 


3C405 


0.056 


1 


27.39 


127 


LMS 


32±3 


36±4 


68±5 


3C330 


0.549 


2 


27.01 


62 


LMS 


82±13 


92±12 


174±18 


3C427.1 


0.572 


2 


27.08 


23 


LMS 


49±5 


67±8 


116±9 


3C55 


0.720 


3 


27.26 


71 


LMS 


73±10 


101±17 


174±20 


3C247 


0.749 


3 


26.89 


16.7 


LPR 


88±13 


97±12 


185±17 


3C265 


0.811 


3 


27.33 


78 


LMS 


82±14 


100±17 


182±22 


3C289 


0.967 


4 


27.28 


11.8 


LPR 


58±6 


53±5 


112±8 


3C268.1 


0.974 


4 


27.46 


46 


LMS 


70±8 


63±7 


133±10 


3C280 


0.996 


4 


27.60 


17.7 


LPR 


52±6 


64±6 


116±8 


3C356 


1.079 


4 


27.43 


75 


LMS 


81±16 


50±10 


131±18 


3C267 


1.144 


4 


27.58 


38 


LMS 


54±9 


51±7 


105±11 


3C68.2 


1.575 


5 


27.85 


22 


LMS 


111±29 


70±10 


181±31 


3C322 


1.681 


6 


27.84 


33 


LMS 


37±4 


47±8 


84±9 


3C239 


1.790 


6 


28.15 


13.0 


LPR 


43±4 


t 


87±6 



V* for only one bridge 
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Table 3. Best Fit f3 min and xl e d Holding < D >/D* Constant, Independent of z. 



case 3. 


n = 1.0, ^ A = o.o 


n a = o.i, 


A = 0.0 


n a = o.i, 


A = 0.9 


Pmin X red 




X red 


Pmin 


X red 


A 


2.60±.35 1.05 


2.30±.30 


0.67 


2.10±.30 


0.69 


B 


2.00±.35 1.84 


2.00±.30 


1.17 


1.95±.30 


1.07 


C 


2.35±.40 0.95 


2.10±.35 


0.67 


1.90±.35 


0.82 


D 


1.70±.30 1.37 


1.75±.30 


1.00 


1.75±.35 


1.11 



a Rcfers to the cases: (A) b = 0.25 without a-z correction, (B) b = 0.25 
with a-z correction, (C) b — 1.0 without a-z correction, (D) b = 1.0 with a-z 
correction. 



-31 - 



Table 4. Best Fit f2 , f3 min) and xled Fitting < D >/D* to a Constant, Independent of z. 







fi A = 




k = 


0, A = 1 - 




case a 


n a 




X red 






X red 


A 
B 
C 
D 


-0.101;" 
-0.351;*° 

o.iol:i8 

-0.05l;« 


2.15±.30 
1.75±.25 
2.10±.35 
1.75±.30 


0.65 
0.95 
0.67 
0.98 


0-201;™ 

0.35l;l 
0.201;^ 


2.25±.30 

n < o 

2.25±.35 
1.80±.35 


0.67 

0.69 
1.07 



a Refers to the cases: (A) b = 0.25 without a-z correction, (B) b = 0.25 
with a-z correction, (C) b — 1.0 without a-z correction, (D) b = 1.0 with 
a-z correction. 
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Table 5. Best Fit (3 min , p(/3 min ), and xled Setting < D >/D, = C(l + 



case a 


n a = 


1.0, n A = 0.0 




Qo = 


0.1, Q A = 0.0 




Qo = 


0.1, n A = 0.9 






p(/3min) 


■^red 




p(/3min) 


Y 2 




p(/^min) 


X 2 j 

^red 


A 


2.25±.35 


-0.75±.30 


0.74 


2.20±.35 


-0.15±.30 


0.71 


2.30±.35 


0.30±.30 


0.70 


B 


1.85±.35 


-0.95±.30 


1.18 


1.85±.35 


-0.40±.30 


1.11 


1.90±.35 


-0.05±.25 


1.15 


C 


2.20±.40 


-0.60±.30 


0.76 


2.15±.40 


0.05±.30 


0.72 


2.25±.40 


0.50±.30 


0.71 


D 


1.75±.35 


-0.60±.30 


1.16 


1.75±.35 


-0.05±.30 


1.08 


1.80±.35 


0.30±.30 


1.11 



a Rcfcrs to the cases: (A) b = 0.25 without a-z correction, (B) b = 0.25 with a-z correction, (C) b = 1.0 without 
a-z correction, (D) b = 1.0 with a-z correction. 
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Fig. I. — Mean lobe-lobe size < D > vs. (1 + z) for high-power, Pns(q = 0) > 3h 2 x 
10 26 W Hz -1 sr" 1 , 3CR radio galaxies for three different cosmologies (see figure legend). 

Fig. 2. — The characteristic length vs. (1 + z) computed without the a-z correction, and 
assuming Q = 0.1, = 0, and b = 0.25 for two choices of (3: (a) (3 = 0.0 and (b) (3 = 2.0. 
Diamonds represent LMS89 sources, and stars represent LPR92 sources. 

Fig. 3.— D/D* vs. (1 + z), assuming Vt = 0.1, fi A = 0, and (3 = 2.0, for b = 0.25 without 
the a-z correction. Since vl is assumed constant over the lifetime of a source, D /D* oc t/t*. 
Diamonds represent LMS89 sources, and stars represent LPR92 sources. The dashed line is 
the best fit to a power-law. 

Fig. 4. — The reduced x 2 as a function of (3 and Q obtained by fitting < D > to a 
constant, independent of redshift, assuming b = 0.25 and no a-z correction: (a) for = 
and (b) fl k = 0. Integer values of the reduced x 2 are labeled in figures; the contours increase 
in steps of 0.25, with the lowest contour at a reduced x 2 of 0.75. Note that the reduced 
X 2 = 1 contour translates to a a of about 2 for both figures, assuming Gaussian errors. 

Fig. 5. — The reduced x 2 as a function of (3 and VL obtained by fitting < D > /D* to a 
constant, independent of redshift, assuming b = 0.25 and with the a-z correction: (a) for 
= and (b) Qk = 0. Integer values of the reduced x 2 are labeled in figures; the contours 
increase in steps of 0.25, with the lowest contour at a reduced x 2 of 1.25. Note that the 
reduced x 2 — 1-5 contour translates to a a of about 2 for both figures, assuming Gaussian 
errors. 
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Fig. 6. — The reduced x 2 as a function of f3 and Q Q obtained by fitting < D > /D t to a 
constant, independent of redshift, assuming b = 1.0 and no a-z correction: (a) for Q A = 
and (b) fl k = 0. Integer values of the reduced x 2 are labeled in figures; the contours increase 
in steps of .25, with the lowest contour at a reduced x 2 of 0.75. Note that the reduced x 2 — 1 
contour translates to a a of about 2 for both figures, assuming Gaussian errors. 



Fig. 7. — The reduced x 2 as a function of (3 and fi G obtained by fitting < D > /D* to a 
constant, independent of redshift, assuming b = 1.0 and with a-z correction: (a) for = 
and (b) fl k = 0. Integer values of the reduced x 2 are labeled in figures; the contours increase 
in steps of 0.25, with the lowest contour at a reduced x 2 °f 1-25. Note that the reduced 
X 2 = 1.5 contour translates to a a of about 2 for both figures, assuming Gaussian errors. 



Fig. 8.— < D > /D* vs. (1 + z), assuming (5 = 2.0, b = 0.25, for three choices of 
cosmology: (a) {Q Q = 1.0, Q A = 0), (b) {Q D = 0.1, Q A = 0), and (c) (Q Q = 0.1, fi A = 0.9). 
Diamonds represent LMS89 sources, and stars represent LPR92 sources. The dashed line 
is the weighted mean in log-space of all < D > /D* in the figure, assuming < D > / D* is 
redshift independent. 
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